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10 \ Abstract . We present an updated overview on the phenomenology of one-loop Higgs 

^^ ■ boson production at Linear Colliders within the general Two-Higgs-Doublet Model 

(2HDM). First we report on the Higgs boson pair production, and associated Higgs-Z 
boson production, at O(ag^) from e'''e~ collisions. These channels furnish cross-sections 

(3 ' i^ the range of 10 — 100 fb for ^/s = 0.5 TeV and exhibit potentially large radiative 

corrections {\5r\ ~ 50%), whose origin can be traced back to the genuine enhancement 
capabilities of the triple Higgs boson self-interactions. Next we consider the loop- 
induced production of a single Higgs boson from direct 77 scattering. We single out 

^ ■ sizable departures from the expected 77 —t- /i rates in the Standard Model, which are 

again correlated to trademark dynamical features of the 2HDM - namely the balance 
of the non-standard Higgs/gauge, Higgs/fermion and Higgs self-interactions leading to 
sizable (destructive) interference effects. This pattern of quantum effects is unmatched 
in the MSSM, and could hence provide distinctive footprints of non-supersymmetric 
Higgs boson physics. Both calculations are revisited within a common, brought-to- 
date framework and include, in particular, the most stringent bounds from unitarity 
and flavor physics. 



1 Introduction and computational framework 

The quest for experimental evidence of the Higgs boson is actively underway at the Tevatron 
and the LHC [T]. Nonetheless, a complete understanding of the Electroweak Symmetry Breaking 
conundrum will not only demand to discover the Higgs boson, but also to precisely measure 
its mass, quantum numbers and interactions to the other particles. A Linear Collider (linac, 
hereafter), such as the ILC or the CLIC, would be the most natural facility to carry this endeavor 
to completion [2] and, perhaps most significantly, to disentangle the Standard Model (SM) Higgs 
mechanism from its many potential extensions. 

The Two-Higgs-Doublet Model (2HDM) [3] constitutes a singularly simple, and yet phenomeno- 
logically very rich example of the latter. The absence of tree-level flavor changing neutral currents 
determines the Higgs/fermion Yukawa couplings and leads to the canonical type-I and type-H 
realizations of the 2HDM [3]. Moreover, it effectively accounts for the low-energy Higgs sec- 
tor of some more fundamental theories, such as the Minimal Super symmetric Standard Model 
(MSSM) [1]. The 2HDM can be fully specified in terms of the masses of the physical Higgs par- 
ticles; the parameter tan^ (the ratio {H2)/{H^) of the two VEV's giving masses to the up- and 
down-like quarks); the mixing angle a between the two CV-even states; and, finally, of one genuine 
Higgs boson self-coupling, usually denoted as A5. We refer the reader to Ref. [5] for full details on 
the model setup, our notation, definitions and various constraints. 

While in the context of the MSSM we expect a panoply of Yukawa, and Yukawa-like, couplings 
of various kinds (including squark interactions with the Higgs bosons), whose phenomenological 
implications have been exploited in the past in a variety of important processes (see e.g. [6]), 
in the case of the general 2HDM we count on alternative mechanisms. Above all, the Higgs 
self-interactions are perhaps the very trademark structure of the 2HDM. Unlike their MSSM coun- 
terparts, the triple (3/i) and quartic (Ah) Higgs self-couplings are not restricted by the gauge 
symmetry, and so they can be potentially enhanced. In favorable circumstances, these enhance- 
ments can translate into highly distinctive signatures of a non-standard, non-supersymmetric Higgs 
sector. Dedicated literature on the topic includes e.g. the tree-level studies on triple Higgs 
boson production, e"'"e~ — )• 3/i [7], inclusive Higgs-pair production through gauge boson fusion, 
e+e" — ^ V*V* — )• 2h + X [8], and the double Higgs-strahlung channels e+e" — ^ hhZ'^ [9]. Also 
significant are the loop-induced single Higgs production, 77 — )• /i [10] and the double Higgs chan- 
nels 77 — )• 2h [H]. They have both been considered in the framework of a photon-photon collider, 
alongside with the complementary radiative decay mode /i — )• 77 |12] . Finally, the impact of these 
3h self-couplings (see e.g. Table H of [5]) has been quantified at the level of radiative corrections 
through the detailed one-loop analysis of the pairwise production of both charged [l3] and neu- 
tral [5] Higgs boson pairs, as well as upon the study of the associated Higgs-strahlung channels 
e^e~ — )• h.^Z^,]i^Z^ at the quantum level |13|ll4j |1|. The effective enhancing power of the Higgs 
self-interactions is subdued, in practice, by a number of experimental constraints and theoretical 
consistency conditions: perturbativity, unitarity and vacuum stability, as well as from the EW 
precision data, the low-energy flavor-physics inputs and the Higgs mass regions ruled out by the 
LEP and Tevatron direct searches - cf. e.g. Ref. |17fl20j . 

2 One-loop Higgs boson production in e''"e~ 

In the following we present a full-fledged one-loop analysis of the CP-conserving pair production 
of neutral Higgs bosons (e'''e~ — t- 2/i = h^A''/H''A'^), alongside with the associated Higgs/Z boson 
(Higgs-strahlung) channels within the 2HDM. On top of the complete set of 0{a^^) corrections. 



^For related work in the context of radiative corrections in Higgs production processes, see e.g. [15]. Phenomeno- 
logical prospects for the LHC have been addressed in |16) . 
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Figure 1: Total cross-section (T(e^e~ — )• h^'A^) (in fb) and relative one-loop correction br (in %) for 
Set A of Higgs boson masses. Left and central panels display these quantities as a function of \fs 
(indicated as Ecm). The results are shown within three different values of A5 at fixed tan/3 = 1.2 
(compatible with the lower tan/5 bound from B^^ ~ B^ data [18]) and for the representative choices 
a = /3 (maximum tree- level coupling) and a = ■njl (fermiophobic limit of h" for type-I 2HDM). In 
the right panel we show their evolution in terms of A5 at y^ = 500 GeV. The shaded areas on the 
left (resp. right) are excluded by unitarity (resp. vacuum stability). 



we also retain the leading 0{a'^^) terms which stem from the (squared of the) 0{X'^i^) Higgs- 
mediated contributions. Renormalization of the SM fields and coupling constants is performed 
in the conventional on-shell scheme in the Feynman gauge [21]. A dedicated extension of the 
on-shell scheme is worked out for the 2HDM Higgs sector in |5]. Phenomenological constraints 
are implemented by interfacing our numerical codes with the packages 2HDMCalc [22], SuperlSO 
[18j and HiggsBounds [23], altogether with several complementary in-house routines. As for the 
algebraic calculation and numerical evaluation of the cross-sections under study, we have employed 
the standard computational software FeynArts, FormCalc and LoopTools [24]. Representative 
choices for the Higgs boson spectrum are sorted out in two sets as indicated in Table [TJ 

Our interest here is basically threefold, namely: i) to seek for regions within the 2HDM pa- 
rameter space sourcing large quantum corrections, which we shall quantify through the ratio 
6r = cr'-^Y'^ ) where a^^' = a — a^^' is the one-loop correction with respect to the tree-level 
value; ii) to evaluate their impact on the overall 2h and hZ production rates (cf. Figs. 1-3); and 
iii) to correlate these effects to the strength of the 3/i self-couplings. 

Worth noticing is that the Higgs/gauge boson couplings {hZZ, hAZ) driving these processes at 
the leading-order are anchored by the gauge symmetry, and hence take the same form in both the 
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Table 1: Sets of 2HDM Higgs boson masses used throughout the calculation. Owing to the flavor constraints on 
Mjj± [IS], Set B is only possible for type-I 2HDM's, whereas Set A is possible for both type-I and type-II. The mass 
sets satisfy the custodial symmetry bound \5p\ < fO~^ - cf. Ref. [5]. 
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Figure 2: Radiative corrections 6r (%) to the total cross-section a{e'^e~ — t- h'^A'^) as a function of 
tan/3 and A5, for Set A of Higgs boson masses, a = /3 and ^/s = 0.5 TeV. The shaded areas in the 
top (resp. bottom) account for the vacuum stability [20] (resp. unitarity |19] ) conditions, while 
the vertical grey band depicts the lower limit tan/3 ~ 1.18 ensuing from B^ — B^ |18j . 



2HDM and the MSSM. Genuine differences between both models should thus be probed through 
the study of quantum effects - among which the enhanced 3h self-interactions of the 2HDM could 
rubber-stamp a very distinctive imprint. 

Let us concentrate on Set A of Higgs boson masses for the analysis of this section. This set is 
possible in both type-I and type-II 2HDM's [3]. (Set B, characterized by lighter Higgs masses, will 
be used for 77 physics in the next section; it is only allowed for type-I models.) The cross-section 
for the particular channel h^A*^ as a function of the center of mass energy ^/s (Ecm), displayed 
for a few values of A5, and also as a function of A5 at fixed ^/s, is shown in Fig. 1, together with 
the relative quantum correction 6r- The rightmost panel of this figure illustrates the expected 
cr ~ o"o + 0{X'l) + ©(As) behavior triggered by the triple Higgs boson self-interactions [5]. The 
cross-sections at one-loop lie in the approximate range of 2 — 15 fb for ^/s = 0.5 TeV - rendering 
10^ — 10^ events per 500 fb~^ of integrated luminosity. The corresponding quantum effects are 
large and positive for ^/s around the nominal startup energy of the ILC, i.e. 0.5 TeV, but become 
rapidly negative for ^/s > 0.6 TeV and stay highly so in the entire 0(1) TeV regime. This sign 
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Figure 3: Total cross-section a{e'^e~ — )■ h^Z'^) (in fb) and relative one-loop correction 6r (in %) 
for Set A of Higgs boson masses. Left and central panels display these quantities as a function of 
y^, for three different values of A5, tan/3 = 1.2 and for the representative choices a = (3 — ir/l and 
a = tt/2. In the right panel we show their evolution in terms of A5 at y/s = 500 GeV. The shaded 
areas on the left (resp. right) are excluded by unitarity (resp. vacuum stability). 



flip in combination with the behavior of the quantum effects on the Higgs-strahlung processes (see 
below) gives an important experimental handle on the physics of Higgs production in the 2HDM. 
Although we have used Set A for the analysis, the behavior of 6r turns out to be fairly independent 
of the details of the Higgs mass spectrum, the particular type of 2HDM and the specific channel 
under analysis (h'^A'^ or H^'A'^). 

A dedicated study of the quantum effects 6r as a function of tan/3 and A5 is explored in Fig. 2, 
at fixed y^ = 500 GeV. We can also appreciate in it the interplay with the unitarity bounds (lower 
area, in yellow) and the vacuum stability conditions (upper area, in purple). Notice that the 
former disallows simultaneously large values of tan/3 and A5, whereas the latter enforces A5 < 1, 
but mainly in the negative range: — 10 < A5 < 0. The largest attainable quantum effects ( 
\6r\ ~ 20 — 60%) are localized in a valley-shaped region centered at tan/3 > 1 deep in the allowed 
A5 < range (cf. Fig. 2). Here a subset of 3h self-couplings becomes substantially augmented - 
their strength growing with ~ IA5I - and stands as a preeminent source of radiative corrections 
via Higgs-boson mediated one-loop corrections to the hAPZ vertex. As a result this interaction 
vertex, which is purely gauge at the tree-level, can be drastically modified at the one-loop order. 
There is, however, the rigid lower bound tan/3 > 1 from B'^ — B^ oscillations [18] (see the vertical 
grey band in Fig. 2), which further restricts the valley-shaped region and finally leaves less than 
half of its original allowance for the largest possible quantum effects. 

In addition to the pairwise production of Higgs bosons we find the more traditional Higgs- 



strahlung channels (e+e~ — )• h.^Z^,I{^Z^) [H], i.e. the 2HDM analog(s) of the so-called Bjorken 
process in the SM [25]. As reported in Fig. 3 (left and central panels), we obtain typical cross 
sections in the ballpark of a(h^Z^) ~ 0(10 — 100) fb, which may undergo substantial (and mostly 
negative) radiative corrections up to order 5r ~ —50% for large 3h self-coupling enhancements 
— these being preferably realized for tan/3 = 0(1) and IA5I ~ 0(10). The trademark negative 
sign of the leading quantum effects can be tracked down to the dominance of the finite wave- 
function corrections to the external Higgs boson fields, this being the only contribution at one loop 
which retains a quadratic dependence on A3/J. The rightmost panel of Fig. 3 nicely illustrates the 
characteristic a '-^ a^ — 0(A5) -|-0(A5) behavior induced by the triple Higgs boson self-interactions, 
which in this case produce dominant (and negative) quantum effects of O^X^) from the wave 
function renormalization as long as the regime IA5I > 1 is well attained |14j . As hinted before, the 
correlation of large negative quantum effects on the Higgs-strahlung channels with the presence of 
significant positive (for ^/s < 500 GeV) or negative (for ^/s > 600 GeV) quantum effects on the 
double Higgs production channels (cf. Fig. 1) could eventually lead to a robust quantum signature 
of (non-supersymmetric) 2HDM physics in a Linear Collider. 

3 One-loop Higgs boson production from 77 

Direct 77 collisions may be realized through Compton backscattering of high energetic laser pulses 
off the original e~^e~ linac beams [26]. This alternative running mode opens up a plethora of 
complementary experimental strategies for a Linear Collider. In particular, it may provide a 
pristine insight into the loop-induced 77/1 coupling - and so to the underlying structure of the 
Higgs sector. This effective interaction ensues from an interplay of gauge boson, fermion and 
charged scalar one-loop contributions |27j . In the 2HDM the charged Higgs- mediated effects, 
which are directly sensitive to the 3/i self-coupling X^^^+h- ? along with the modified Higgs/fermion 
and Higgs/gauge boson interactions, are responsible for a highly characteristic phenomenological 
pattern. To illustrate it, we shall concentrate on the following quantities: i) the total (unpolarized 
and averaged) cross-section {a^^^h){s), which results from the convolution of the "hard" scattering 
cross-section (7(77 — t- h) with the photon luminosity distribution (describing the effective e — )■ 7 
conversion of the primary linac beam). For the latter we use the parametrization included in 
the standard package CompAZ [28]); ii) the relative strength r of the effective 77/1 interaction 
normalized to the SM (upon identifying Mh^^^^ with M^a), namely r = g-y-yh/g^^n^ ^ly which we 
aim at quantifying the departure of the genuine dynamical features of the 2HDM with respect to 
the SM. 

The results of our numerical analysis of r are displayed in Fig. [H A first observation is that, 
unsurprisingly, the profile of g^^\^o is mainly modulated by the strength of the 3h coupling Ajjojj+jj- . 
This is a consequence of a destructive interference between the contributions from the loop diagrams 
triggered by the Higgs boson self-interaction h'^H^H" and those induced on the one hand by the 
gauge bosons and those driven by the Yukawa interactions with neutral Higgs bosons and fermions 
on the other (cf. Fig. 2 of |10j). The impact of such interference is well visible for wide areas 
of the parameter space where the strong departures from r ~ 1 can be traded to suppressions of 
order 40 — 60% for the effective g^^h interaction as compared to the SM. Away from these depleted 
domains, maximum cross-sections of order a ~ 0(10) fb - viz. up to a few thousand events per 
500fb~^ - are attainable for both h*^ and H^ (cf. Fig. [5] for the dependence of the cross section with 
sin a and tan/3). Interestingly enough, the optimal production rates are nicely complementary for 
the CV-even channels 77 — >■ h'' and 77 — )• H*^, as a result of the inverse correlation of the respective 
A^jj+jj- self-couplings (see Table II of [5]) - and hence of the dominant interference effects. In 
contrast, owing to the CP-odd nature of A", the 77 A*^ channel appears to be rather featureless and 
has a milder numerical impact. 
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Figure 4: Effective 77h'^ coupling strength in the 2HDM normalized to the SM, r = g^-yh" / 9j^h ' 
in terms of sin a and tan/3. The yellow strips signal the lower and upper bounds stemming from 
unitarity, while the grey vertical band ensues from the B'^ — B^ constraints at the 3o" level. 



Besides the 3/i self-couplings, further 2HDM mechanisms could contribute, at least in principle, 
to shift the r ratio from its canonical value r = 1. For example, r > 1 could be achieved for 
tan/3 ~ 0.2 — 0.3 (cf. Fig. U]) as a result of the enhanced Higgs-top (~ l/tan/3) Yukawa coupling. 
In practice, however, the B^ — B^ bounds exclude this possibility. Significant departures from the 
SM may arise as well within large A5 scenarios, such as those analyzed in [10]. This possibility 
is nonetheless highly disfavored if the unitarity conditions are included in their most restrictive 
version [19] . 

All in all, trademark hints of a 2HDM structure may emerge from 77 — )■ h, mainly through 
a missing number of events with respect to the SM predictions - as long as the overall rates are 
still large enough to be efficiently discriminated from the dominant background process, 77 — )• hh. 
We have shown that this situation can still be realized in sizeable regions of the parameter space, 
provided the 3/i self-couplings take on moderate values (viz. \^h ~ O(IO^) GeV) preserving 
unitarity. In contrast, in the MSSM the genuine supersymmetric (slepton/squark- mediated) con- 
tributions to g-y-yh can only induce rather tempered quantum effects as compared to the general 
2HDM, the main reason being the absence of potentially large 2>h self-couplings, and hence the 
lack of a mechanism able to prompt the characteristic interference pattern that we have identified 
above. Although alternative enhancing effects on 77 — )• /i within the MSSM are possible [29], e.g. 
through (light) stop- mediated loops with large trilinear {At) couplings and sizable mass splittings 
between their chiral components ti,t2, they are nevertheless comparably weaker. In fact, these 
effects are always pulled down by inverse powers of the SUSY breaking scale and are further limited 
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by the stringent limits on the squark and Higgs boson masses, as well as from i?-meson physics 
- cf. the recent Ref. [3D] for a fully updated MSSM analysis and a comparison with the general 
2HDM case. The bottom-line is that mild deviations from r = 1 of order —5% characterize the 
typical MSSM scenarios. In this sense, it is worth emphasizing that, even if a pattern of the sort 
r < 1 would overlap with the 2HDM predictions for small 2>h self-couplings (viz. A3/J ~ 10 GeV), 
the correlation of 77 — t- h*^ and 77 — t- H*^ could still help to disentangle both models. Indeed, in 
the 2HDM both CT^-even channels could simultaneously yield a ~ 1 — 10 fb (see Fig. 5), whilst 
such situation is definitely precluded in the MSSM owing to the SUSY restrictions on the Higgs 
boson mass splittings - see [30] for details. 



4 Discussion and Conclusions 

In this work, we have described several phenomenogical aspects of one-loop Higgs boson production 
processes within the general 2HDM at the future Linear Colliders. We have revisited our previous 
results and cast them into a common, fully updated framework, including the most recent set of 
theoretical and experimental constraints presently available in the literature - most significantly 



those stemming from unitarity and flavor physics. Our brought-to-date analysis keeps on high- 
lighting the truly instrumental role reserved to the future linac facilities, and also the great degree 
of complementarity of the e+e~ and 77 running modes. 

We have provided detailed, quantum-corrected predictions for the exclusive pairwise production 
of Higgs bosons e^e~ — )■ 2/i as well as for the Higgs-strahlung channels e+e~ — )• hZ at 0{a^^) and 
leading 0{a^^). In the case of 2h production, we have shown that the radiative corrections can 
reach the level of \5r\ ~ 50% for enhanced 2>h self-couplings - mostly for low tan/? and IA5I ~ 0(10) 
- and hence give rise to a substantial positive boost with respect to the tree-level expectations 
(around ^/s ~ 500 GeV) or suppression (for -y/i > 600 GeV). In the case of the hZ final states the 
quantum effects can be of the same order of magnitude, but they are negative for essentially the 
whole y/s range. Optimal rates for these processes lie in the ballpark of 0(10 — 100) fb and can 
be attained for typical Higgs masses in the 100 — 300 GeV range. 

No less crucial is the role of the 3/i self-couplings in the production of a single Higgs boson from 
direct 77 collisions. A trademark suppression of (7(77 — )• h) with respect to the SM predictions 
is singled out, and its origin traced back to the interference effects between the different one-loop 
contributions, these being critically modulated by the strength of the 3/i self-couplings /iH~''H~. In 
regions of the parameter space for which this depletion is moderate (viz. A^^jj+jj- ~ 0(10 — 100) 
GeV), still significant cross-sections up to a ~ 10 fb may be retrieved for both h*^ and H*^. Prom 
the experimental point of view, the opportunities for accessing this kind of final states are deemed 
to be excellent. The decay signatures of the neutral CP-even states would essentially boil down to 
h*^ —7- bb/T'^T~ or h*^ — t- VV — t- 4/, 2l+$j>: ^^ of them allowing for a comfortable tagging in the clean 
linac environment. Interestingly enough, the described phenomenology is particularly distinctive 
of a non-supersymmetric 2HDM structure. In the MSSM, potential enhancements cannot be 
triggered by the Higgs self-interactions - which are anchored by the gauge couplings - but instead 
by the Yukawa interactions of the Higgs bosons and the sfermions. These give rise, in general, 
to rather tempered quantum effects, as compared to the sizable corrections that are spotlighted 
for the 2HDM. Likewise, the conditions that SUSY dictates on the Higgs spectrum may be of 
relevance here. Por instance, genuine indication of non-standard, non-SUSY Higgs physics may 
come from the simultaneous observation of 77 — )• h" and 77 — )• H*^; both channels may yield O(IO^) 
events per 500fb~^ - a situation which could never be ascribed to the MSSM, as the mass splitting 
between the two Higgs bosons is enforced to be much larger. Pinally, the combined analysis of 
these signatures together with complementary multi-Higgs production channels (cf. Ref. [Sj) could 
unveil a characteristic pattern of signatures for different values of ^/s. If confirmed, it would point 
to a non-standard, non-supersymmetric origin. Conversely, if the two CP-even states would be 
produced at measurable rates differing by, say, one order of magnitude, this could be compatible 
with MSSM Higgs physics (see [30] for details), but it would require a detailed dijet invariant mass 
reconstruction to resolve the spectrum and check if it is compatible with the MSSM constraints. 
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